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Abstract - From a scale analysis of hydrodynamic phenomena having a significant action on the drift 
of an object in coastal ocean waters, we deduce equations modeling the associated hydrodynamic fields over 
a time period of several weeks. These models are essentially non linear hyperbolic systems of PDE involving 
a small parameter. Then from the models we extract a simplified and nevertheless typical one for which we 
prove that its classical solution exists on a time interval which is independent of the small parameter. We 
then show that the solution weak—* converges as the small parameter goes to zero and we characterize the 
equation satisfied by the weak—* limit. 
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1 Introduction 

This paper is part of a work program concerning the modeling of object drift in near coastal ocean 
waters over a several week time period. 

The final target of this program is to develop methods to forecast the drift of things like con- 
tainers, lost objects or oil spill over long periods of time in near coastal ocean areas. Such methods 
would be of interest for services in charge of maritime safety, environmental studies or pollution 
impact assessment. To reach this target, several research topics need to be further investigated. 
For instance, improvements are needed in the field of the numerical methods to simulate long term 
drift, in the modeling and simulation of the near coastal ocean waters, in the understanding of 
ocean-object and ocean-spill interactions and, of course, in the integration of all those aspects to 
move toward a coherent theory. 

In the previous paper of this work program, Ailliot, Frenod and Monbet [2| considered the 
numerical method facet. We built a numerical method coupling a two scale expansion method, 
explored in [Hj, and a stochastic wind simulator, in the spirit of [TJ [31 [Mli in order to estimate 
probability of events that may happen to the considered object, such as running aground in a given 
area. In [5], the simplified model was supposed to describe the object dynamic in the ocean. It 
involved an ocean velocity field which was decomposed into a sum of a velocity due to the tide wave 
and of a perturbation. Both of them were periodic of the tide period and with modulated amplitude 
and the fields used for the numerical simulations were not realistic. 



The present paper deals with the modeling of near coastal ocean. The sea velocity and the 
fluctuation of the sea level due to the tide wave are well known in many coastal areas of the world. 
Those main fields are perturbed by fields with a smaller order of magnitude having a net long-term 
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result. Such perturbations, which are produced by meteorological factors, propagate and interact 
with the main fields. The precise aim of the present paper is to make a first step toward the set 
up of a modeling procedure in order to establish partial differential equation systems describing the 
evolution of those perturbations and to suggest ideas to solve them. 

The paper is organized as follows. In section [21 we summarize the main mathematical results. 

Then, in section [3l we set up the previously evoked modeling procedure. The models which 
are built in this section are deduced from the Shallow Water Equations via a scale analysis of the 
geophysical phenomena concerned and of the geometrical size of the concerned domain. As for 
the size of the coastal domain and the characteristic order of magnitude of the wind velocity, we 
consider several possibilities giving rise to several models. They are all essentially hyperbolic systems 
of partial differential equations with a singular perturbation involving a small parameter. 

From those models, in section |4l we extract a simplified one which is nevertheless typical. For 
it, adapting classical methods for hyperbolic systems (see Kato [T^, Majda [55], Klainerman and 
Majda [HI [19], Schochet [3T1[32], and Metivier and Schochet [23]), we prove that its classical solution 
exits with a time existence independent of the small parameter. 

In section [5] using a homogenisation method (see Tartar |35j , Bensoussan, Lions and Papanico- 
laou [5], Sanchez-Palencia [29] and Lions [20]), we set that this classical solution weak—* converges 
to a function. Keeping within the mind frame of Frenod [S], Frenod and Hamdache [TU], Frenod, 
Raviart and Sonnendriicker [TT], Joly, Metivier and Ranch [TCj or Schochet [33] we finally look for 
the form of this function and establish the equations allowing for its computation. 

Finally, in section [6] we conclude and give some perspectives. 

acknowledgments - The authors want to thank Joanna Roppers for proofreading the manuscript. 



2 Results 

In this section we present the main results. We first present one of the models, involving a small 
parameter, set out in this paper. Then we state a Theorem claiming the existence of the solution 
to a simplified version of this model. Finally, we exhibit the asymptotic behavior of this solution as 
the small parameter goes to zero. 

The model we present now, and which is set out among others in section|3l describes the evolution, 
over a several months time period, of the perturbation of the sea velocity and of the sea level in 
an ocean domain above the continental shelf at a latitude about 45° and with stormy weather 
conditions. 

The small parameter involved in this model is the ratio tide duration on observation time scale. 
The first one is about 13 hours and the second is about three months. Hence the involved small 
parameter is: 

e=—. (2.1) 
200 ^ ' 

Variables and fields involved are all rescaled; rescaled meaning that the order of magnitude of those 
variables and fields is one and that they have no physical dimension. 

The rescaled velocity of the sea M and the water depth H induced by the tide wave are considered 
as known and periodic with modulated amplitude. In other words, t being the rescaled time and x 
the rescaled position, 

M(t,x) = M(i, -,x) and H{t,x.) = Hit, -,x), (2.2) 

£ £ 

where M and TL are regular functions and where 9 ^ (M(i, 6, x),H{t, 6, x)) is 1— periodic. 
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The model says that the total sea velocity, expressed in km/h, writes 0.5(M + eN) and that the 
total sea level is -^{E + 2e'H + 2e^T) , where E is the rescaled mean sea level and where N and X 
are rescaled perturbations. Moreover, (T, N) is solution to 

^ + V{^E + 2n)-N+ i^E + 2n)V ■ N + 2(Vi) • M + 2X(V • M) 3) 

+2e{(yi) ■ N + i(V • N)) = 0, 

— + 2(VN)M + 2(VM)N + 2£(VN)N + — + — VX - 13e*AM - ISe^AN 

_^g^4(VM)VGB + 2£H) _ (VM)vi _ ^^^,(yNME + 2en) 

E + 2eH + 2e^i E + 2e'H + 2e^i E + 2en + 2e^i 

-26.^^XN)V^ 

E + 2eH + 2e2j [2^) 
„ 1 1 



E+2sn+2e^I M + 3 E+2en+2e^I ji;^ _ 



s 1 + !^{E + 2eH + 2e^I) 1 + ^ {E + 2sH + 2e'^I) 
1 , 1 1 , 

6 E+2eH+2e^I - M) - 6e E+2eH+2e^I jyj 

1 + ^{E + 2671 + 26^1) e 1 + ^{E + 2eH + 2e^i) 

In this system, W is the rescaled wind velocity, N-'- = N2, Ni), A stands for the laplacian, V- for 
the divergence operator and V stands for the gradient of scalar fields and for the Jacobian matrix 
of bi-dimensional fields. 

Motivated by this system, we consider a simplified version of it which consists in considering that 
the ocean bottom is flat, i.e. E = 1, in forgetting all the power of e greater than 1 and in setting all 
constants to 1: 

f)i 1 

— + (V7Y) • N + (- + 7i:)(V • N) + (Vi) • M + X(V • M) + e((VZ) • N + Z(V • N)) = 0, (2-5) 
at e 

^ + (VN)M + (VM)N + £(VN)N + -N^ + -VI = W. (2.6) 

ot e £ 

In this system, t £ [0, T], x = {xi,X2) G K^. The unknowns are 1 = l{t, x) and N = N(t, x). Their 
evolution is influenced by M and H for which we assume (12.21) and by W for which we also assume 

W(t,x)=W(t,J,x), (2.7) 

with function W regular and with 9 i—>- W(i, 0,x) being 1— periodic. This assumption is not really 
convenient for real wind velocity time series but is comfortable from a mathematical point of view 
(see Ailliot, Frenod and Monbet [5] for a more detailed discussion). Moreover, we equip this system 
with the following initial conditions 

ilt=o=io, N|t=o = No, (2.8) 
and we can claim the following Theorem. 

Theorem 2.1 Under assumptions / fg.g|) and |g.7| j, 1/(10, No) £ {H^{R'^))^ with s > 3, then there 
exists a time T, not depending on e, such that the classical solution (I, N) £ (^([0, T], (iJ*(R2))^) n 
(Ci([0,r],(i7"-i(M2))3)) of and exits and is unique. Moreover this solution sat- 

isfies 

sup ||(X,N)||, <c, (2.9) 

te[o,T] 
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for a constant c not depending on e, where \\ \\s stands for the norm in {H'^{M?))'^ . 

Concerning the asymptotic behavior of (Z, N), as e goes to zero, we have the fohowing result. 

Theorem 2.2 Under the assumptions of Theorem \2.1[ there exist functions X = X{t,x) e C([0,T], 
i/*(R2)) and N = N(i,x) e C([0, T], (iJ"(]R2))2)^ gy^^h that as e goes to 0, the solution (X, N) of 
^K^, figi and weak-* converges to (J, N) in L°°([0,T], (i?^(]R2))3). Moreover, I and N 

are linked by 

dl dl 
Ni(i,x) = - — (i,x), N2(<,x) = — (i,x), (2.10) 
0x2 0x1 

and X is solution to 



■M-VI 



dt dxi dxi 



8x18x2) V 8x1 



8x2 8x2 

'8M2 81 \ .^/8M^ 81 \ ^/8M^ 81 \ ^/8M^ 81 



i-P^) 

V 8x1 8x2 ' V 8x2 8x1 / V dxi 8x2 / V 



8x2 8x1 



8x1 8x1 8x2 8x2 

dWt 8W^ 



8x2 8x1 

where M = / M.d6, H — I TLdO and W = / Wd6', and equipped with initial conditions 



, (2.11) 



3 Models 

In this section, we first consider a reference modeL It consists in removing the ocean level and the 
ocean velocity which are induced by the tide wave from the Shallow Water Equations. This gives 
rise to a system of equations governing the time evolution of the ocean level perturbation and of 
the ocean velocity perturbation. Then, we analyse the scale of the variables and fields involved in 
the problem we want to describe. Rescaling the reference model in view of this scale analysis finally 
yields the desired models. 



3.1 Reference model 

It is generally admitted that the evolution of the ocean level h = h{t, x) and of the ocean velocity 
m = m(t, x) is well described by the following Shallow Water Equations 

8h 

— + V{h-hb)-m.+ {h-hb)V-m. = 0, (3.1) 

s , \ „, . (Vm) V(/i - /lb) 7;-^ 

— + (Vm)m + /m-L + gVh - cAm - c> + ILJi^---m = 

8t ^ ' h-hb l + -(h-hb) 

^ ' (3.2) 

.^H;!'' , , (W~m) + F, 
1 + 7(/i - hb) 
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Figure 1: Fields h, hb, E, H and X 



equipped with ad-hoc initial and boundary conditions. This system was introduced by Saint- Venant 
[27] . For an exhaustive explanation concerning ocean modeling and the construction of this model we 
refer for instance to Pedlosky pS', Nihoul [25 , Lions, Temam and Wang pT, Stoker Whitham 
[37] or Johnson [15] . For a deduction of the Shallow Water Model taking into account viscosity, being 
able to model the consequences of wind and bottom friction actions, which is considered here, we refer 
to Gerbeau and Perthame [14]. In system ()3.ip - ()3.2|) . hi, = /ib(x) is the depth of the ocean bottom, 
/ is the Coriolis parameter, g is the gravity acceleration and c is the water viscosity. The friction 
coefficient on the bottom is k and the air-water friction coefficient is /i. Lastly, W = W(t, x) is the 
wind velocity and F may take into account the action of other meteorological factors like atmospheric 
pressure. 

Now we isolate the action of the tide wave. In other words, we consider that the ocean depth 
variation 7i = Ti{t, x) around the mean water height E = E{x.) and the ocean velocity M = M(t, x) 
which are induced by the tide wave are known. We consider that {E + ?{, M) is the solution to 

^ V(£: -f 7Y) • M + (£: + 7i:)V • M = 0, (3.3) 
^ + {VM)M + fM^+gV{E + n + hb) = 0, (3.4) 

with initial and boundary conditions imposed by the tide wave. 

A brief parameter size analysis, that will be confirmed in the next sections, shows that the terms 

— cAm — c— — r— ^ — + . , ^7,'''' , . m have a very small influence on the sea movement. Hence we 

have chosen to put them in the equations for the perturbations hereafter. 

Now we introduce the perturbations I and N which are defined such that h = hi, + E + H+ T 
and m = M + N (see Figured]). 

Replacing h and m by these expressions in p.ip - p.2p and removing the terms appearing in 
p.3p - p.4p leads to the equations for T and N. Then we obtain the following reference system which 
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is the starting point of our scale analysis. 

di .... .... 

— + V{E + ■ N + + 7^) (V • N) + (VZ) ■ M + J(V • M) + (VX) • N + X(V • N) = 0, (3.5) 



'dt 



(VN)M + (VM)N + (VN)N + /N^ + .gVi - cAM - cAN 

(VM)V(i; + 7i:) (VM)VZ (VN)V(i; + 7i) (VN)VI 
—c- — = c = — c c- 



E+n+i E+n+i E+n+i E+n+i 



E+H+I jyj ^ E + n+I jyj- 



(3.6) 



1 + ^^{E + H + T) l + '^iE + H+I) 
E+n+i ixkr ivrN E+n+i 



:^(W - M) 



i + f(£; + w + x) i + f(i; + ?^ + x) 

3.2 Rescaled variables and fields 

We introduce a reference time Z, two reference lengths L and 7. Those reference values, as well as 
the other ones introduced hereafter, will represent characteristic values (mean or maximum values 
for example) of the physical quantities under consideration. We consider the rescaled variables t' 
and x' =(x']^, X2) expressing time and position in unit t, L and I. They are defined as 

t — tt' , xi — Lx'i and X2 = 1^2- (3-7) 

If the reference values are chosen as evoked above, the order of magnitude of the rescaled variables are 
1. Then we define M and TV as the characteristic velocity of the tide wave and its perturbation; E the 
characteristic value of the mean water depth, H the characteristic tidal range and / the characteristic 
value of its perturbation. W is the characteristic wind velocity and F the characteristic scale of the 
field F. The rescaled fields have the following definitions: 

M'(i',x') = ^M(Zt',Lx;,Ix^), N'(i',x') = ^N(tt',Ix;,7x^), (3.8) 

i?'(x') = ii?(Lx'i,7x'2), (3.9) 

7i;'(t',x') = =7T:(tt',rx'i,Ix'2), i'{t\x') = ji{tt',Lx[~ix'2), (3.10) 

W'(i',x') = =^W(tt',Ix;,Ix^), F'(i',x') = 4F(tt',Ixl,7x^). (3.11) 
W F 

Lastly, we introduce ZU the tide wave frequency and we assume that M and Ti, are l/ZZJ— periodic 
function with modulated amplitude. In other words, we assume 

M'(t',x') = M'(t',cJZi',x'), 7Y'(i',x') = n'{t',ojtt',x'), (3.12) 

where 9 ^ {M'{t',e,x'),n'it',e,x')) is 1-periodic. 

From system (|3.5p - p.6l) . we deduce the following rescaled equations for N', and T' with known 
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fields M',E\ W and F': 



dX' H Nt 



( E dE' 

L E dE' 



dW 

dx 





V I dx' ^ dx'o 7 



+ !'{ 



dx'^ 



dx'n 



m 
T 



_dx'^ 
L dl' 



Mt 



dx'^ 
LdT 



■M' 



(3.13) 



(9N' 



Mt 

T 



/^N' LdN' 



/dM' LdM' 

T 54" 



N' 



NtrdN' LdN' 



L V ' I dx'2 



- N' + ftN 



at I 
+= = 
NL 



dx'^ 
LdT 
\ld7j 



ct M d^M' 



9x2^ 9a;'i^ dx'2 



72) 



/dm' LdM' 
V dx'i ' I 9a;2 



ct M 



( dE^ '^dfe_ \ 
_dx{ Wdx{ 
L dE' Hdn' 

Kl^d^^^^d^'J 



Jf N 



E' + &H' + Li' 



/dM' LdM' 
V dx'i ' I dx'2 



ct M I 



_dx'^ 
L dX' 
\ldZj 



/dN' LdN'\ 

y'd^'J'd^) 



ct 



I dE' Hdn' \ 
L dE' Hdn' 

KVd^^^^d^'J 



L^ NE 



E' 



En' 



E' 



En'- 



'dN' LdN'^ 
dx'i ' I dx'2 ' 



ct I 



( dX^\ 
_dx'^ 
LdX' 



l'e 



E' 



H 



n' 



KtM 

+== 



E N 1 + !iE{E' + En' + LX') 



—M' + = 



Ei + nE^E' + En' + LX') 



-N' 



E I + iiE(^E' + Eii' + Lf') N N E i + bE^^e' + En' + Lf') 



Ft 

N' + =F'. 

N 



(3.14) 



3.3 Parameter size and rescaled equations 

In this subsection, we fix the characteristic values. As set out in the previous subsection, we shall 
choose mean values or maximal values of the concerned physical quantities. The parameter t is the 
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observation time scale. We consider that it is about several months. Then we set 

i - 100 days - 2400/i, (3.15) 
beside this ZJ is the tide frequency, meaning l/uJ is the tide duration, i.e.: 

L ~ 13/i. (3.16) 

Hence we exhibit a small parameter: 

1 1 
^ 200' 

Then, we make a strong assumption, which is that N and X are really perturbations 
words, we consider that 

N I 
= ^ = ^e. 

Concerning the Coriolis parameter, for latitudes about 45° / ~ Ti/day ^ 4 10^^/s, then ft ^ 
7r/2e. Concerning the other parameters of physical meaning, several choices are possible, according 
to the turbulence action, the nature of the ocean bottom, the shape of the ocean free surface 
and so on. We focus on one of those choices, being aware that others, that would lead to other 
models, are also reasonable. For the viscosity, we chose the value of the water viscosity at 20° C, 
i.e. c ~ 10~^cm^/s ~ IQ"'' km? / day , then ci ^ 10~^krr?. Concerning the friction coefficients, the 
bottom friction coefficient is k ~ 10~'*m/s ~ lO^^km/day and the air-water friction coefficient is 
H ~ 10~^rn/s ~ 10~'^km/day. Those values are consistent with the ones used, for instance, in 
Dawson and Proft [7] or Gerbeau and Perthame [TJ. Then kI ~ 1km, fit ^ 10~^km, k/c ~ 10^ /km 
and fi/c ~ 10^ /km. We also have g ^ 10m/ ^ lO^km/day^ and gt ^ lO^km/day. 



We now turn to the ratios determining the asymptotic analysis we have to realize. Having in 
mind our final target, i.e. the drift of things in the ocean over long time periods, we notice that 
such a drift may take place relatively far from the coast, above the continental shelf. It may also 
take place in a large and relatively closed bay, with a long residence time of the ocean water. Such a 
domain will be called coastal zone. As was the case in 1999 / 2000 for the Erika oil slick along the 
French Atlantic coast, the drift may occur for weeks along a thin layer following the coast. Those 
remarks guide the choices concerning the geometrical assumptions we consider further. 

L and I represent the characteristic lengths of the domain where the drift takes place and M /U 
the characteristic distance the water covers in the tide duration. Following Salomon and Breton 
[28] . Cetina, Rajar and Povinec [36], Bao, Gao and Yan [4J or Cai, Huang and Long |6j, we can 
state that this distance is about a few kilometers (from 5 to 20) in the cases we are interested in. If 
the domain under consideration is a continental shelf, the characteristic sea water velocity M is 
about 0.5km/h and then, we have 

M 

— - 5km, (3.19) 

Ll) 

and if we set L ^ 500fcm, E ~ 300to and H ^ 3m, then 

^ - 2e, = - 2e. (3.20) 
L E 

We also have iV ^ sM ^ 2.5 IQ^^km/h ~ 6 10~^km/day, then we get gt/7^ ^ 1.7 10^. Since 
7 ~ e!ff ~ 1.5 lO^^m, we obtain I/L ~ 3 10"*. Hence 

oil 1 , , 

= =-50-—. (3.21) 
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Moreover 

ct 10" 



-2 25 104 



13£^ (3.22) 



Kt 1 kE . 0.8 

= ^-^3.3^3, -^J5 10^^- (3.23) 

d^l^^3.3 10-^6e, ^^3 10^^1:^. (3.24) 

Concerning the wind velocity, when the weather is calm, 10km /h is a relevant characteristic value, 
while 100km /h may be a good choice in stormy conditions. Hence, we shall consider 

^ 0.5 1 

W ^ 10 ^ 20' 

in calm weather regime, and, 

E^^^e, (3.26) 
100 ' ^ ^ 

in storm regime. Expressing now the following ratios 

m - M- m JliN 

^^tuj^, ^-^ = , 3.27) 
L L L L M 

and moreover setting Ft ^ N and removing the ' for clarity, we can write the rescaled equation 
(|3.13[) - (|3.14|) in the case of a continental shelf: 

5 + Vi-E + 2H)) ■N+{-E + 2n)V ■ N + 2(VJ) • M + 2i(V • M) 

at £ £ {6.10) 

+2e((vi) • N + i(V • N)) = 0, 

— + 2(VN)M + 2(VM)N + 2e(VN)N + — N-^ + — VX - 13e*AM - 13e^AN 
ot 2e 4e 

_^g^4(VM)V(^+2£H) _ (VM)VJ _ ^^^^ (VN)VGB + 2£H) 

S + 2eH + 2e2i + 2e7i: + 2e2x + 2e7T: + 2£2i 

^^'^ E + 2£'H + 2£^J (3-29) 
1 1 



S+2e-H+2£2l M + 3 E+2en+2e^I p;;^ _ 



£ 1+ ^j§.{E + 2£H + 2£^I) l+^jiiE + 2£n + 2£^I) 

1 1 



E+2en+2e^I - M) - 6£ E+2en+2e^X _N + F. 



l + ^f{E + 2£H + 2£^I) 1 + ^f{E + 2£n + 2£^I) 

where 7 = 20 = l/(10e) in calm weather regime. In storm regime, which is what it is supposed for 
model (|2.3p - (|2.4p presented in the introduction, 7 = 1/e. 

If the domain is a coastal zone, M 1km/ h 

M 

— ~ 10km, (3.30) 

UJ 

and we set L ^ 5km^ E ^ 50m et H ^ 10m. In this case 

¥ H 1 , , 

^-2, =--. 3.31 
L E 5 ^ ^ 
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We also have iV ~ eM ~ 1.2 IQ-^km/day, then we get gt/N ~ 8 10^. Since J ~ eiJ ~ 5 IQ-^m 
we obtain I/L ~ 10~^. Hence 



Moreover 



Kt 







~ 8 10^ 


0.2 








NL 










d 


10-5 , 






1 

0.05 


~ 20 ~ 


1 

TOe' 


C 


5 10^ 


1 

10£2 




10-2 
0.05 


0.2, 


c 


50 ~ 


1 

4e' 



Concerning the wind velocity, we have 

M 1 

W " To' 

in calm weather regime, and, in storm regime 



(3.32) 

(3.33) 
(3.34) 
(3.35) 

(3.36) 



= - - 2e. (3.37) 
100 ^ ^ 

Hence, the rescaled equation reads in this case: 

dT 2 ~ ~ 2 ~ ~ 2 ~ 

— + -{y{5E + n))-N+ -{5E + H)V • N + -(VJ) • M + -J(V • M) 

at e £ e £ (6.6a) 

+2{Vi) •N + 2i(V-N) =0, 

^ + -(VN)M + -(VM)N + 2(VN)N + — N-^ + - O.Se^AM - O.Se^AN 

at £ £ 2e £-^ 

_Q (VM)VCB+^ _ ^ (VM)VT _ ^ (VN)VCB+^ 

-0.1.3 (VN)Vi 



(3.39) 



0.2 ^ ^+^^+§5 _(^W - -M) - 0.2 , ^+i^+t^ _N + F. 

l + ^(^+iW+|J)'2£ e i + i^+feJ) 

where 7/2 = 10 = l/20£ in still weather and 7/2 = l/2£ in stormy weather. 

We will give the name of coastal layer to a domain having the following characteristics 

M 

— ~ 10A:to, (3.40) 

CO 

and L ~ 500km, I ~ 5kni, E ~ 50rn et H 10m and then 

M 7 II 1 

^ ~ 4£, = ~ 2£, = ~ -. (3.41) 
L L E 5 ^ ' 
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In this case, we have N ~ eM ~ 1.2 10 ^km/day, then we get gt/N ~ 8 10*. Since / ~ sH 
5 10~^m, we obtain I/L ~ 10~^. Hence 



9t I 0.4 

= = ~ 80 . 

N L s 



Moreover 



ct 



10- 



1 



Kt 



E 0.05 



20. 



25 104 
1 kE 



13e^ 



Ht IQ- 



E 0.05 



lOe' 



0.2, 



5 W 



c 
c 



10£2' 



50 . 



4e 



(3.42) 

(3.43) 
(3.44) 
(3.45) 



The considerations concerning the wind velocity are the same as in the case of a coastal zone. The 
the rescaled equation for a coastal layer writes: 



dl 
'dt 



oxi ox I 

2 aE an 

\e 8x2 9x2 / 



N + (20£' + W)( 



9Ni 1 9N2 



dx\ 2e 8x2 



( ^ \ 

dxi 

1 di 



M 



dxi 

\ dX2J 



(3.46) 



9a; 1 



dxo 



0, 
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dt \ dxi ' e dx2 / V dxi ' e 8x2 / V dxi ' 8x2 J 2e 



,dM 2dM 



+ - 



0.4 



/ 9J \ 

dxi 
1 di 



(13e' 



dxx^ 



4 dx2^ ' 



tdm 1 9M\ 

\9a;i ' 2e 




/5M 1 aM\ 

V dx\ ' 2e 



) 



1 d± 
\2^9^/ 



/9N 1 9N\ 

V 9xi ' 2e 9x2 ' 



US'" 




/9N 1 9N 

V9xi ' 2£ 9x2^ 




+ 



lOe 1 



10; 



0.2- 



_^(^W- -M) -0.2- , 



— = 



(3.47) 



4 Existence 

4.1 Simplified system for continental shelf 

In this section, we focus on one of the models introduced in the previous section, and we explore 
some of its mathematical properties. 

More precisely, we consider a simplified version of system p.28p - (|3.29p which consists in consid- 
ering that the ocean bottom is flat, i.e. E = 1, in forgetting all the power of e greater than 1 and 
in setting all constants to 1. Then we obtain system (|2.5p - (|2.6p and we prove an existence result for 
the solution of this system. 

Although the results given in sections [4] and [5] are specific to the model (|2.5p - ()2.6|) . we expect that 
similar methods could be used to prove similar results for the other models introduced in section [3l 
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4.2 Proof of Theorem [211] 

Setting u = (i,N) = (X, Ni, N2), u^^ = (0,N^) and introducing 

/Mi + eNi -^U + el \ 

1 

Mi+eNi 

p 

V Mi+eNi/ 

+ eNs -+'H + ei\ 



B^Ct, -,x, eu) 
e 



B^{t, - ,X, £U) 

e 



V 



M2 + eNa 

M2 + / 



and 



i^(t,-,x, u) 

e 



equations (|2.5p - (|2.6p read 



And, introducing 



axi oxg^ OTi 0x2 

Wi - - — Ni + - — N2 

ax\ 0x2 

W2 - {— — Ni + — — N2 

axi 0x2 



du du , 3u 1 

ot 0x1 0x2 e 



£ 







A^{t, -,x, eu) 
e 



V 

/ Ml + eNi 

1 + eTY + en 




o\ 

1 
Ij 



Ml + eNi 
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Ml 
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M2 


+ eN2 




M2 + 






X, eu) 




1 


+ eTY + e^f 



eN2 








I 












M2 
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0^ 






^0 


l\ 


S^ = 


1 
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0/ 
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equation (j4.4p yields the following symmetric hyperbolic system: 

,ndu .T du du 1 1 9u 1 9u 1 , ^ , „ ^ , , „n 

^1^ + + + -S^— + -S^— + -u^ = = A^'F. 4.9 

Hence applying Kato [17] or Majda [S^, we deduce that for any e, under the assumptions of Theorem 
12. li the classical solution of (|2.5[) , (|2.6p and (|2.8p exists and is unique on a time interval, what remains 
to prove is that this time interval does not depend on e. 

Before going further in the proof, we may observe the following differences between equation 
(|4.9p and the type of problems, also depending on a small parameter, studied in Klainerman and 
Majda [18l|T9], Schochet [311 [30l [32] , and Metivier and Schochet ^3]. Some of those differences 
simplify the problem: the nonlinearity in A'^ are functions of only e^I and in A^ and A"^ of only eu. 
Some others make the results proved by those authors unable to be applied directly: A'^, A} and A? 
depend on i/e and the singular term u^/e involves the function u itself and not order 1 derivatives 
of it. Nonetheless, the now classical calculus procedures carried out in the concerned papers may be 
followed in order to obtain the right estimates allowing for the conclusion. We sketch the concerned 
computations hereafter. 

We set a — (q;i,q;2) G with lal — a\ ^ a-i < s and D"u — — — — — --. Applying D" to 

0x1 0x2 

equation (|4.9p yields 



with 



A"—— + A^— + A^— + -S^— + -S^— + -GD"u^ = F'^, 4.10 

at 0x1 0x2 e 0x1 e 0x2 e 



= D^Fo - [D'^,A^-]n-[D^,A'—]n~[D",A'—]n, (4.11) 
ot 0x1 0x2 

] standing for the classical commutator. 

Multiplying equation (|4.10p by 2D"u, , integrating on and noticing that 

dD"u d( [ A'D^u-D'^udx) o^ 

AO^!^ . D'^u dx = - / ^i?"u . D-udx, (4.12) 

dt dt J dt ^ \ J 

dD^Vi ( f diAW^u- D^u) \ f d(A^) 
A^— D^udx^ / -5^ ^dxl- / -V^D"u-D"udx 



'^^'^^''D"u-D"udx (4.13) 



dxj 



S3^^^±-li.D"udx=-2 S^-- i?"dx = 0, (4.14) 

OXj J OXj 



for j = 1, 2, and 

we obtain 

A°D"u- D^udy: 



2 J (£i"u)-L • D^'udx = 0, (4.15) 



dt I dt 



f !M^D"u- D''udx+ [ ^li^L»°u-L>"udx+ [F^-D^udx. (4.16) 
J dxi J dx2 J 



14 



For all the estimates to come, all the constants which are needed are called c. Since the dependency 
of with respect to t/e is done through eH, and since s > 3, we can deduce that for any t and x, 



dt 



dt 



at 



(4.17) 



where || ||s-i stands for the norm in H'^ ^(R^). The time derivatives of the other entries of are 
zero. Hence the first term of the right hand side of (|4.16p may be estimated 



dt 



ax, 



<c{l + e'\\^l^,)\\D'^nro<c(.l + s^\\j^l_,)\\u 



dt 



(4.18) 



where || ||s stands for the norm in (_ff*(R^))-^ and || ||o for the norm in (_L^(]R^))^. Concerning the 
entries of A* for i = 1, 2, 



diAl^it,^,^,eu)) 



dxj 



< c(l + £ 



. dl 
' 9a;,; 



) < c(l + e sup 



dxi 
du 



+ e sup 



dJ 

dxi 



<c(l + £|| — 1|^_^) <c(l + e||u||J. (4.19) 



Hence 



^^D'^u-D^ud^ 

dxj 



<c(l+£|lu|U)||i?"u||^<c(l+£|lu|U)||u|| 



(4.20) 



The last term of ()4.16p is left to estimate . For this, we first notice that D"Fq is the sum of 
controlled coefficients multiplied by D^u (possibly multiplied by or or . . . ) with f3 = (/3i, /?2) G 
such that (3 < a (i.e. /3i < ai and (52 < a2)- Hence 



< c(l + ^ P^llo) p"u||o < c(l + ||u|U) ||u|| 



(4.21) 



Secondly [D", A^g|^]u is the sum of controlled coefficients multiplied by D^u (possibly multiplied 
by e or £^ or ... ) and themselves multiplied by D'^u with < a, 7 < a and /3 + 7 < a + (1, 0) 
which implies |/3| + I7I < |a| + 1. When < s - 1 and I7I < s - 1 since \(3\ + I7I + 1 < |a| + 2 < 2s 
we deduce that D^n ■ D'^u G L'^{M?) with ||_D'^u • Z?'''u||o < cl|it||^ by a classical calculus inequality 
that can be for instance found in the Appendix of Schochet [30]. When \a\ = s, \j3\ = s and I7I = 1 
we have sup^g^s l^'^'u] < ||u||s. Then D^'u • D'^'u e L'^{R^). When laj = s, \(3\ = 1 and I7I = s 
suPxGR2 l^'^ul < l|u||s and then Dl^u ■ D^u e L^{M.'^). As the same can be done with [D°',A^^]u, 
we deduce 



\D'^,A^—]u-\D",A^—]u] •L>"udx 
0x1 0x1 ' 



<c(l + ||u||. + |ju||^)P"u|| 
<c(l + ||u|U + ||i 



(4.22) 



Finally, [CjA^^Ju is a sum of controlled coefficients multiplied by eD^u (possibly multiplied by 
e or £^ or . . . ) and themselves multiplied by D''^ with 7 < a and (0,0) < f3 < a, that implies 
\P\ > and + < \a\. When < s-1 and I7I < s-2, since |/3| + |7| + 1 < |a| + l < s+(s-l), 
applying classical calculus inequalities, we obtain D^u ■ D'^^ & L^(R^), with ||Z?'^u • D'^^Wq < 

chliyitll^^' < When \a\ = s, \P\ = s and I7I = we have 
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suPxGM^ ID'^^I < Then D^u-D^^ G L'^{^^). When |a| = s, = 1 and I7I = s - 1, we 

get sup^gR2 \Dl^\i\ < \\u\\s and D^u ■ ^ G L^{R^). Hence, we deduce 



< c(l +£ 



du 



dt 



Us U 



(4.23) 



Using inequahties (|4.17p - (I4.23P and summing (|4.16p for a < s, we obtain 



E 

\a\<s 



dt 



<gi(||u|U,e||^|U_i), 



(4.24) 



for a function gi not depending on s. 

Derivating system (|4.10p with respect to t, we get 

ot 0x1 0x2 e 0x1 e 0x2 e ot 

= ^^D-^~^D-^~^D-^ + ^. (4.25) 
dt dt dt dxi dt dx2 dt 

Every previously established estimate remains valid. Moreover, we can set that the entries A].i of 
A* for i = 1, 2 satisfy 



d(4;(t,|,X,£u)) 



dt 



.1 II II ||9U|| , 

<^(- + ^ll"IL + ^ll5^IL-i). 



(4.26) 



and using, when it is necessary, the same classical calculus procedure as above, we can show that 

dt 



is in i2(K2) ^i^-jj 



dF° 



dt 



,1 

<c - + 1 
e 



II l|2 l|l9U|| l|(?U||2 ,|| I 

-ll"IL + ^II^IL-i + ^ll5^IL-i)ll"l 



(4.27) 



Then, multiplying equation (|4.25p by 2e_D"(e^) and following the previous method, we obtain 

</.OZ,«(.f,-Z,.(.f,.,) 



\a\<s-l 



dt 



<g2(I|u|U,e|||^|U„i), 



(4.28) 



for a function §2 not depending on e. 

As a conclusion, estimates (|4.24p - (|4.28p together with the fact that 



{J2 J A°D"u- D'^udx 



1/2 



(4.29) 



\a\<.s 

is a norm equivalent to || ||s, allows us to set that 





du 




d(\Ms+s 


dt 







dt 



<5(l|u|| 



du 



dt 



.) 



(4.30) 



for a function g not depending on e and then that the time interval on which the classical solution 
of (|2.5p . (|2.6p and (|2.8p exits does not depend on e. 



Finally, estimate (|2.9p is a direct consequence of (|4.30p . This ends the proof of Theorem 12. II 
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5 Asymptotic behavior: proof of Theorem 12.21 

In order to deduce the asymptotic behavior as e goes to of (Z, N) we use the method, developed 
in Tartar [5S], Frenod [8j and Frenod and Hamdache [TU] and used in Frenod and Sonnendriicker 
]T2[ [T5] , which consists in setting a weak formulation of system (|2.5p - (|2.6I) or of its equivalent 
form (|4.10p . Passing then to the limit using the weak—* convergence allows us to set a constraint 
equation. This constraint equation imposes a form to (X, N). Using test functions satisfying the 
constraint equation allows us finally to deduce system ()2.1ip . 

Estimate ((^ yields the weak-* convergence of u = (I,N) to U = (I,N) in L°°{[0,T], 
H^{M.^)^), up to a subsequence, as e goes to 0. 

Multiplying symmetric hyperbolic system (|4.9p by test functions 5'(t,x) being regular, R'^ — 
valued, and with compact support in [0, T) x R^, and integrating yields 

u • — — + + + -J. + -S^^ + -S^^ + -4-^ ]dtd^ 

\ ot e oO oxi OX2 £ oxi e 0x2 e J 

^[ [ A"F ■■iidtdx+ [ uo • A°*(0,0,-)dx. (5.1) 
Multiplying (|5.2p by e and passing to the limit yields, since A^ converges to / and OA'-' /d9 to 0, 

-T [ ru.(s^^-,S^^ + ^-de)dtd.^O, (5.2) 
Jo Jk^Jo V "xi dx2 J 



which is the weak formulation of 



+ + (5.3) 

0x1 0x2 



V • N = 0, + VI = 0. (5.4) 

From this constraint equation, we deduce the form of {X, N) given by (|2.10p . 

For any regular function cp we define the test function 5" satisfying the constraint equation by 

vl'i(t,x) =^(t,x) vl/2(t,x) --|^(t,x) M/3(f,x) = |^(t,x). (5.5) 

0x2 0x1 

Using this function in (|5.ip cancels terms containing 1/e factors. Since A° converges to /, 1 /e dA^ /dO 
weak—* converges to 0, A^ weak—* converges to Midd I and A^ weak—* converges to Jp 'M.2dd I 
passing to the limit yields 







f ^ f8^^^iJ,^M,do)^^^iJ,^M2de)^^ 

Jr2 \ Ot dxi 8x2 ) 

= [ [ [ Fde--^dtdy:+ [ Uq • *(0, O^x. (5.6) 
Jo Jr2 Jo Jffi2 
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or, using expressions of U and F, 



: V dt dxi 0X2 



8x2 V dt dxi 8x2 

dj d(f^Mide)^ d(f^M2de)S^^ 

oxi V OT OTi aa;2 / 

^ djj^nde) dx ^ dij^Hde) dx ^ , d{j^M,de) ^ dij^'M2d9) 

V dxi 8x2 dx2 dxi dxi 8x2 / 

\ Jq 0x1 0x2 0x2 oxi J 0x2 

vjo i^^;! 0x2 0x2 0x1 / 0x1 

+ f io^-(No)i-^ + (No)2|^dx. (5.7) 

We have here a weak formulation of (|2.1ip . Since this equation is Hnear, it is easy to show that its 
solution is unique. From this, we can finally deduce that the whole sequence u weak—* converges 
to U as £ ^ 0, ending the proof. 



6 Conclusion and perspectives 

In this paper, we set out equations modeling the long term evolution of the perturbation X of the 
ocean free surface elevation and of the perturbation N of the velocity field. Because of the tide 
wave, those models contain and generate oscillations with high frequency. If numerical simulations 
of near coastal ocean waters are needed, directly using those models could be very expensive because 
of the oscillations. Nevertheless, the result given in Theorem 12.21 . which says that (T, N) weak—* 
converges to (Z, N), suggests a way to use those models for numerical simulations of near coastal 
ocean waters. As an intuitive interpretation of it we could say that 

X{t,x.) is close to X{t,x.) and N(t, x) is close to N(t,x). (6.1) 

Hence, since equations (j2.1ip - (|2.12p neither contain nor generate oscillations with frequency 1/e, 
we can solve it using a numerical method involving a time step which does not need to be small 
compared with e. Hence solving (|2.1ip - (|2.12p and then reconstructing (X, N) via (|2.10p in place of 
solving (|2.5p - (|2.6p - (|2.8p directly may give good results in a shorter computational time. 

Among the tasks listed at the beginning of this paper, we realized significant steps in the direction 
of building numerical methods in our previous paper [2] and in coastal ocean modeling over long 
time periods in the present one. 

The next step, we shall provide in a forthcoming paper, will consist in using the models set out 
here in order to compute the ocean fields in a real near coastal ocean area and to couple this to 
the numerical method proposed in [2'. This will allow us to make forecasts in a real coastal ocean 
region. 
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